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Abstract 

We present a theoretical and experimental study of spin precession in the presence of both a 
static and an orthogonal oscillating magnetic field, which is non- resonant, not harmonically related 
to the Larmor precession and of arbitrary strength. Due to the intrinsic non-linearity of the system, 
previous models that account only for the simple sinusoidal case cannot be applied. We suggest 
an alternative approach and develop a model that closely agrees with experimental data produced 
by an optical-pumping atomic magnetometer. We demonstrate that an appropriately designed 
non-harmonic field makes it possible to extract a linear response to a weak dc transverse field, 
despite the scalar nature of the magnetometer, which normally causes a much weaker, second-order 
^ \ response. 

PACS numbers: 32.30.Dx, 32.10.Dk, 32.80.Xx, 
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I. INTRODUCTION 



The precession of spins in a static homogeneous field is a well known problem of classic as 
well as quantum physics. A second widely studied system is derived from the previous one 
when an oscillating field perpendicular to the static one is added. This is the configuration 
used in typical magnetic resonance (MR) experiments, and for this reason, it has been 
extensively studied. In MR setups, however, the transverse field oscillates at a frequency 
near-resonant to the Larmor frequency set by the static field, and has usually an amplitude 
much smaller than the static one. The resonant condition greatly facilitates the task of 
modeling the system: the oscillating field is seen as a superposition of two counter-rotating 
(and usually weak) fields, one of which appears to be (quasi-) static in a reference frame 
which rotates at the Larmor frequency around the static field, while the other has negligible 
effects. This is well-known in textbooks as the rotating wave approximation (RWA). 

In this work we study the system described above, but in conditions which make the 
RWA unsuitable. In particular, we consider the cases in which the transverse field has a 
general periodic time dependence containing Fourier components at frequencies much larger 
than the Larmor frequency. In fact, we re-examine, extend, and generalize the results 
reported decades ago by S.Haroche and co-workers where similarly an analytical model 
has been developed and applied to interpret experimental data from optical detection of 
atomic precession. 

The results presented here are modeled with the Larmor equation and discussed using a 
perturbative approach based on the Magnus expansion of the time-development operator. 
Compared to previous works on this subject, our approach makes it easy to model the effects 
of a transverse field oscillating with an arbitrary waveform. 

It seems worth stressing that the subject has implications for the research area of dynamic 
localization J2, 3)], which concerns the effect of a time-periodic field on quantum systems such 



as charged particles in crystals |4| and Bose Einstein condensates 

The experimental implementation is feasible, provided that a slow precession occurs 
around a weak static field, for instance with ultra-low-field NMR setups and in atomic 
spin precession experiments. This latter is the nature of the setup used in our case, and the 
reported experimental results concern measurements performed with an optical-pumping 
atomic magnetometer operated at a /zT static field range. The corresponding Larmor fre- 
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quency is in the kHz range, making it easy to apply a homogeneous transverse field oscillating 
at much higher frequencies, even with a much larger amplitude. 

Our apparatus detects the time-dependent Faraday rotation of the polarization plane 
of a weak probe laser beam, which is directly related to a component of the macroscopic 
magnetization of atoms. The analyzed magnetization component is the one parallel to the 
direction of the oscillating field, which coincides with the probe beam propagation axis. In 
the presence of a strong transverse oscillating field the three components of the magnetization 
evolve in time with a complex behaviour, which is far from a simple precession. Nevertheless 
the measured component maintains its similarity to the simple precession behaviour. In 
other terms, the motion around the direction (z) of the static field is deeply affected by 
a strong non-resonantly oscillating field along x, which mainly modifies the dynamics of 
the magnetization components M y and M z , but leaves M x evolving (approximately) with 
a simple harmonic law. This quasi-harmonic evolution of M x makes it possible to define 
a modified (dressed) Larmor frequency even when the oscillating field is so strong that no 
precession motion could be recognized by simply observing the time evolution of the whole 
M vector. 

The paper is organized as follows: in Section |TT] we describe the experimental setup used, 
in Section II III we discuss the theoretical model, and Section IIVI contains the discussion and 
the conclusions. 

II. EXPERIMENTAL SETUP 

The basic structure of the experimental setup is that of an optical atomic magnetometer 
(see Fig. [TJ). Using a resonant circularly polarized laser light, alkali atoms are pumped 
in a specific Zeeman sublevel and thus the medium is magnetized along the wave-vector 
(x) direction. The atomic medium is placed in a homogeneous magnetic field {z direction) 
transverse to the pump light propagation. Consequently the Zeeman sublevel where atoms 
are pumped evolves in time, giving rise to the magnetization precession in the xy plane, 
around the static magnetic field. 

This precession is detected by analyzing the Faraday rotation of a weak, linearly po- 
larized, probe beam which propagates in the xy plane (in our case it is collinear to the 
pump beam). The linear probe beam polarization can be decomposed in two a + and a~ 
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counter-rotating circularly polarized components, which experience variation of the respec- 
tive refractive indexes due to the atomic spin precession. These two circular components 
recombine with a time-dependent relative dephasing in a linear polarization rotated by an 
angle oscillating at the Larmor frequency. Beside the different dispersions, the <j + and o~ 
components suffer also different absorptions, but this aspect is made of little relevance by 
tuning the probe laser away from the resonance center. 

To contrast the relaxation processes that cause a progressive decay of the induced mag- 
netization, and hence of the detected signal, the pump light is applied periodically, syn- 
chronously with the precession of the atomic spins, i.e. resonantly with the Larmor frequency 
set by the static magnetic field. 

In our setup, the pump radiation is generated by a single-mode, distributed feedback, 
pigtailed DFB diode laser, whose optical frequency is controlled through the junction current 
and is periodically made resonant (or non resonant) to the Dl transition of Cs, at 894 nm. 
It illuminates the atomic sample at an intensity of about 0.5 mW/cm 2 . The probe radiation 
is generated by a single-mode, Fabry-Perot, pig-tailed diode laser resonant with the D2 line 
of Cs, at 852 nm. It illuminates the sample at much weaker intensity, being attenuated at 
a level of about half /iW/cm 2 , in order to negligibly perturb the sample. 

Polarization maintaining (PM) fibres bring the radiation of the two lasers to a coupler 
(via an on-fibre attenuator in the case of the probe) which mixes them in the PM fibre illu- 
minating the sample. In proximity of the atomic sample, the two radiations emerging from 
the fibre end are collimated by a lens, their linear polarization is reinforced by a polarizing 
cube. Then, a phase-plate acts differently on the polarization of both the radiations, render- 
ing the pump radiation circular, while leaving the probe one linear. This double effect is 
achieved thanks to a multi-order wave plate made of a 480/im quartz layer, which introduces 
|f| for the crossed linear polarization components of the two wavelengths a (5+1/4)Adi and 
a 5A^>2 relative delays, respectively. Such method of tailored multi-order waveplate method 
is cuite similar to that applied in the magnetometric setup (based on Rb vapors) described 
in [7|. 

After interacting with the atomic sample, the pump radiation is blocked by an interfer- 
ential filter and the probe beam polarization is analyzed by a balanced polarimeter made 
of a 45° oriented Rochon and a couple of photodetectors whose photocurrent difference is 
amplified by a low-noise transimpedance amplifier. 
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Figure 1: Schematics of the experimental setup. The probe laser is tuned to the D2-Cs line 
@852nm and unmodulated; the pump laser is tuned to the Dl-Cs line @894nm and is modulated 
synchronously with the atomic precession. Both are mixed and coupled to the sensor though 
polarization maintaining fibres (pmf). The probe is attenuated down to fiW level before the Cs 
cell, the pump is stopped by an interferential filter (Int fltr) after the cell. A multi-order waveplate 
renders the pump polarization circular, while leaving the probe polarization linear. The Cs cell 
contains 23 Torr of N2 as a buffer gas, which quenches the fluorescence and slows down the Cs 
atom motion, making it diffusive. The cell is slightly heated to increase Cs vapour density. The 
transmitted probe laser polarization is analyzed by a balanced polarimeter made of a Rochon 
polarizer oriented at 45° with respect to the probe polarization axis and two photodiodes. A 
differential transimpedance amplifier converts the photocurrent to a voltage signal, which is then 
digitally acquired. That signal reproduces the Faraday rotation, associated to the time-dependent 
orientation of the precessing atomic spins. 
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The polarimetric signal is either analyzed as a function of the laser-modulation frequency 
to characterize the atomic magnetic resonance, or used to generate a signal for driving the 
laser optical frequency, thus closing a loop which makes the system a self-oscillator. In the 
first case, the precession frequency is estimated with a best fit procedure, while in the second 
case it is evaluated as the self-oscillation frequency, provided that the loop delay is properly 
compensated. In practice, the closed-loop, self-oscillating operation helps in shortening the 
measurement time, but may introduce systematic errors due to inappropriate choice of the 
loop delay. On the contrary, operating in the open loop/scanned frequency/best fit procedure 
requires a longer time but improves the accuracy in the Larmor frequency determination. 
All the plots reported in Sec. IIVI were obtained in this more accurate open-loop operation 
mode. With the used static field intensities, the precession frequency ranges from hundreds 
of Hz to several kHz. 

The atomic sample is made of a high quality gas-buffered cell containing a droplet of Cs 
and 23 Torr of N2. The buffer gas is suitable for slowing down the thermal motion of Cs 
atoms, making it diffusive and also quenches any stray fluorescence light. Following several 
optimization campaigns devoted to improve the sensitivity of the setup in magnetometric 
applications, the pump laser is made resonant (at the right time) with the | 1 5'i/2, F g = 4) — > 
\ l P\/2j hyperfine components of the Dl line, while the analyzed ground state is the F g = 4 
and the detection is made tuning the probe laser about 2 GHz on the blue wing of the line 
made of the triplet \ 1 S 1/2 ,F g = 4) -» | 2 P 3 /2)- The pump radiation is broadly modulated with 
a slightly asymmetric signal, so that, besides the synchronous excitation of the | 1 S'1/2 , F g = 4) 
level, it also spends some time in resonance with the \ 1 Si/2,F g — 3} — >■ | 2 -P3/2) triplet. 

The Cs cell is inserted in an anti-inductive heating coil by which the temperature is 
increased up to about 50°C in order to reach a good compromise for having a large signal 
(an optimal optical depth) without introducing an excessive broadening due to the increase 
of spin-exchange collision rate. A second coaxial, inductive coil makes it possible to apply 
a homogeneous time dependent transverse field which is thus parallel to the optical axis. 
This second coil is supplied by a digital waveform generator which provides an arbitrary 
periodic voltage signal. This signal is designed to produce current (and hence x-oriented 
magnetic field) Fourier components with appropriate phases and amplitudes, in accordance 
to the measured complex impedance of the coil. 

The whole experiment is conducted in a magnetically clean volume, where a set of three 
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mutually perpendicular, large size Helmholtz coils and a set of quadrupoles make it possible 
to control the three components of the magnetic field and reduce their inhomogeneities. In 
particular, the horizontal components (x, y) of the Earth field are fully compensated while 
the vertical one (z), partially compensated down to a few /xT, determines the static field. 

The system is here used for an absolute frequency measurement. Nevertheless, the ap- 
paratus is designed to perform differential measurements as well. In fact, the fibre coupler 
is a 2 x 2 50% device, meaning that the two radiations are mixed in equal ratios and made 
available in two output PM fibres. Thus, the apparatus can easily be adapted to work with 
a dual sensor, suitable for gradiometry and for differential magnetic measurements. The ac- 
curacy in the absolute frequency measurement (order of 1 Hz) is largely in excess compared 
to the needs of this work, and is limited by the environmental magnetic noise. 

III. MODEL 

The aim of this section is to describe and characterize the precession of spins in the 
presence of a time dependent magnetic field orthogonally oriented with respect to the static 
one, using a special form of the perturbation theory for the time evolution operator. 

As stated above, taking the z axis in the direction of the static field and the x axis along 
the alternating one, i.e. B = [B x <f)(u)t), 0, B z ), the Larmor equations for the magnetization 
are 

m x =-n z M y 

< M y = tt z M x - Q x <f>(ut) M z (1) 

]\d Z = Q X (f)(Ujt) My 

Here Q xz = jB X}Z (7 is the gyromagnetic ratio) and <p(ut) is a real and periodic function 
with a well-behaved Fourier expansion. Notice that hereafter the magnetization is treated 
classically, but the same results hold true if the Mi, i = x,y,z are treated as non-commuting 
quantum mechanical operators. 
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We are interested in the case of ( = Q z /u <C 1, so it is natural to use perturbation theory. 
To keep track of the perturbation order, we introduce the parameter e and let £ — > e( in the 
intermediate formulas, we set then e = 1 at the end of the calculations. We can thus rewrite 
eq. d2J) as 



—M=[e(A z + £<J ) (T)A x ] M, 



(3) 



where £ = Q x /u. Notice that introducing the commutator A y = [A z , A x ], the three matrices 
A&, k = x,y, z constitute essentially the three-dimensional real representation of the angular 
momentum quantum operators. 

The magnetization M(r) can be obtained acting with the time evolution operator U(r) 
as M(t) = C/(t)M(0) on the initial value M(0). In the spirit of perturbation theory, it is 
useful to factorize U in the "interaction" representation as 



/ 



U (r) = exp 



{r')dr'A 



UAt) 



1 
cos $ — sin $ 
sin $ cos $ 



UAt) 



(4) 



which implicitly defines Uj{t). Moreover we set 

$(r) = £ f (f)(T')dr'. 



(5) 
Jo 

Let's comment on the structure of eq.(|4]): the first factor represents a rotation of an angle $ 
around the first axis. This rotation relates the laboratory frame to the rotating frame, thus 
$ is the Larmor angle for the precession around the oscillating field. In the rotating frame, 
the evolution operator satisfies 



dUi 
dr 



e ( e -Hr)A xAze Hr)A x Vi = £ £ [ cog $A Z + sin $ Ay] XJi = e(A^ Uj 



whose solution, following the Magnus expansion device |8 



can be written as 



The exponent satisfies the non-linear equation 8| 
dW 

. a A" ' - 

2 



< 4 7) - ~ {W, eC 4 7) ] + 4 i W > \- W > < 4 J) ] ] + - 



(6) 



(7) 



(8) 



dr ~^" z 2 l " 1 ^" z J ' 12 
which can be solved perturbatively. In fact expanding W as W = Xm>i e ™ ^™ one S e ^ s ^ ne 
first order contribution 



W 1 {t) = C Re(*(r)) A 2 + C Im(*(r)) A/ 



(9) 



where we introduced the quantity 



(10) 



After some algebra one finds 

( cos(C|*|) -cos/3 sin(C|^|) sin/3 sin(C|*|) 

Ui= cos/3 sin(C|*|) cos 2 /3cos(C|^|) + sin 2 /3 sin(2/3) sin 2 (C|^|/2) 

v -sin/3 sin(C|*|) sin(2/3) sin 2 (C|^|/2) sin 2 /3 cos(C|*|) + cos 2 /3 J 

which represents a rotation of an angle (\^\ around the axis 



n = (o, sin ft, cos /3 j • 



(12) 



Considering the initial magnetization prepared along the x axis M(0) = (M x (0), 0, 0) and 
using the eqs. (ITT]) and (@J one finds 



M x (t) = M x {0) cos(C|*(r)|), 



(13) 



clarifying that \l/ is the Larmor angle around the axis n. Deriving this quantity, the first- 
order correct frequency of the signal experimentally analyzed can be evaluated 



a L ) = ^ (CI*(t)I) 



d 



(14) 



Notice that, on general ground, if one needs only the argument of the cosine in f TT3"j) . it can 
be more easily obtained by calculating the eigenvalues of W\. 
Another case of practical interest is the following 



d 



—M = [e ( A z + e 2 £ A x + £ 0(r) A x ] M, 



(15) 



which models the case of a very small (residual) static magnetic field superimposed to the 
alternating one. Precisely we assume that the dimensionless quantities £, ( and £o satisfy 



£ = lB^/u « 0(1) C = £lB z /u « 1 £ = e^B^/u « C, 



(16) 



showing the need of second-order perturbation theory to study the effect of £o- Repeating 
the above derivation we get the same W\ while the second-order contribution W2 becomes 



W 2 (r) 



A x = X (r) A, 



(17) 
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The geometric meaning of \ is that in the second order approximation the rotation described 
in eq. (jH]) does not occur around n, but around a direction having also a small x component. 
Finding the eigenvalues of eW\ + e 2 W2 makes it possible to determine the argument of the 
cosine 

e L (r) = e^C 2 \^(r)\ 2 + e^x(r) 2 (18) 



and thus eq. (I14p generalizes to 



n?) = jblW (19) 



A. Determination of \& and x 

We are interested in periodic signals, so that 4>{r), $(t) and exp(i$) have a well behaved 
Fourier expansion. More precisely we consider signals <fi{r) with a zero mean so that 

oo oo 

e i*(T) = G n e inr = \G n \e^ nT+e "\ (20) 

n=— oo n=— oo 

then it is easy to obtain 

[ T e iHT') dr > = G r + V — (e inT -l) =G t + 9 (t). (21) 

where the implicitly defined function g(r) is a limited and oscillating function. When r ^> 1 
(or t ^> 1/w) we can neglect g(r) obtaining ^(t) = Got, and thus we find a dressed Larmor 
frequency 

n^ = n t \G \. (22) 

To access the second-order corrections we need the integral 

f ^(r')e-^dr' = \G \ 2 r 2 /2 + n{a - b) + / (r) + r/i(r) (23) 
Jo 

where /o(t) an d /i(t) are periodic and limited functions whose explicit form is not important 
here. The parameters a and b are defined as 



|2 



aeV^ b = y\^L. (24) 

Neglecting the oscillating terms in eq. (123|) as before and taking the imaginary part we get 

X(r) =t[£ + C 2 (6 - He(o)) /2] = r + C 2 c/2) = r Xo . (25) 
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The dressed Larmor frequency eq. fU9p can thus be evaluated approximating the square 
root in eq. (jTBJ) as 

G °\ + ^TTT [(^o/C) 2 + ^oc + C 2 c 2 /4] , (26) 

where now we can consider the last terms as corrections to the first-order result (122]) . Notice 
how the coefficients G n , n^O are encoded in the c parameter. 



n 



a 



\G(] 



1 + 



2C 2 |G | 



B. Pure sinusoidal signal 

This case has already been studied and reported in literature. S. Haroche et al., in 
a pioneer paper had already demonstrated the possibility of solving the simple case of a 
harmonic transverse field [1] with a different approach based on quantum mechanics. In 
that study, the authors quantized the oscillating field as a bosonic mode (rf photons) and 
used standard first-order perturbation theory to evaluate the energy shifts. Then, in the 
(semi-classical) limit of a large boson number, they obtained the same result reported here 
below. That method would run into difficulties for generic time-periodic signals. On the 
other hand, let us examine how quickly the results are obtained using our approach. Starting 
: rom 0(r) = cos r we easily get $(r) = £ sin r and the coefficients G n are the Bessel functions 

G n = \J n (0\e ie " (27) 

being the phases 9 n either or tt according to the sign of J n . We thus re-obtain the Haroche 
et al. result 

n L = n z \Jo(0\- (28) 

Notice that in case of a general harmonic signal 0(r) = cos(r — r ) one has 

$(r) =e[sin(r-r ) +sinr ] (29) 

and thus 

G n = e ^ sinT0 \J n (01 e^""™ 7 " ) (30) 
giving for |Go| the same result as before. 
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C. Non-harmonic signals 



Having found the general expression (I2"2"j) for the dressed Larmor frequency, we can study 
more complicated cases. Let's first introduce a higher harmonic component in the signal, 
indicating with r the ratio between its amplitude and the fundamental one i.e. 



which gives 



[t) = cos r + r cos(kr + <p) 



$(r) = £ sin r + — (sin(A;r + ip) — sin tp) 



Repeating the same steps as before, one gets 



G, 



-i ($r/k) simp 



^2 J n-k P (0 Jpiirjk) e 



(31) 



(32) 



(33) 



p=— oo 



which is a generalized Bessel function [10]. Thus the leading contribution (eq. ([22]) ) of the 
Larmor frequency becomes 



a 



F(£,r,k,<p) = 



J -k P (0 Jp(£r/k) 



p=—oo 



(34) 



Another highly non-harmonic signal which gives analytical results is the square wave 



$(r)=£ 



< r < 2na 
a(2n — r)/(l — a) 2na < r < 2n 

In fact in this case the integration outlined in eq. (12"0"|) can be done trivially obtaining 

r _ -ferfq sin(7rea) 
'jo — e , 

and thus the dressed Larmor frequency results 

Ql |sin7r^a| 



(35) 



(36) 



(37) 



For general periodic signals the G n coefficients can be obtained analytically as multiple 
sums of products of ordinary Bessel functions. However, this approach soon becomes cum- 
bersome and a numerical evaluation of the G n , using for instance the Fast Fourier Transform, 
is a viable alternative. 
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IV. DISCUSSION 



The main issue of the paper is to study the effect of an alternating, non-resonant magnetic 
field oriented perpendicularly to the static one causing a reduction of the precession rate of 
spins. Such a reduction has a non-monotone dependence on the amplitude of the alternating 
field, and, under specific conditions, the precession frequency can take a zero value. This 
phenomenon has been previously observed and theoretically interpreted in the case of a 
harmonic signal 

The model introduced in this work is suitable to the description of the simple case of a 
harmonically oscillating field but also makes it possible to easily extend the analysis to the 
non-harmonic case. 

The comparison of theoretical predictions and experimental results is made taking into 
account several operating conditions and varying different parameters. We preliminarily 
test the consistency with known results obtained in the simple case of a transverse field 
oscillating with a harmonic law. Then we test the strongly non-harmonic case of a square- 
wave field. A more detailed analysis is presented for the case of a two-Fourier-component 
field, where we analyse the effect as a function of the amplitudes and of the relative phase of 
the two components. We also devote attention to the cases when the frequency ratio is even 
(k = 2) or odd (k — 3), pointing out the symmetries discussed in the appendix [A] Finally 
we demonstrate the existence of a linear response to a weak dc field, which has implications 
in the design of atomic magnetometers with vectorial response. 

We show in Fig. [2] the comparison of experimental data obtained with the apparatus 
described in Sec. [Ill and the curve (namely |Jq|) given by the model in the case of a harmonic 
field. In this figure (as in the following ones) the amplitude of the oscillating field is expressed 
in terms of the dimensionless quantity £ defined (see also Sec. IIHI eq. (J3])) by £ = £l x /u> = 
jBx/u. The frequency reduction effect is also given in natural dimensionless units, reporting 
the ratio between the dressed and the unperturbed precession frequencies. It is worth 
noticing that if the same curves are plotted in terms absolute quantities (i.e. plotting 
the dressed precession frequency as a function of B x ), passing from a spinning species to 
another having a larger gyromagnetic factor, the vertical axis scales up while the horizontal 
one scales down. This makes it possible to identify conditions for matching the dressed 
precession frequencies of different species, which is of great interest for possible applications, 
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see for example 



H] and references therein. 
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Figure 2: Ratio between the dressed and unperturbed precession frequency as a function of the 
amplitude of the oscillating field (expressed in the dimensionless quantity £ = jB x /uj) in the case 
of a harmonic transverse field, according to eq. (|28p the plot reproduces the behaviour of | Jo (£)|. 

As a first example of a non-harmonic modulation, we consider the case of a symmetric 
(50% duty cycle) square wave. In this case the model predicts a behaviour given by a sine 
modulus (eq. (135]) ). The prediction is in good agreement with the experimental observation, 
as appears evident in Fig. [3j 

Square wave signal 
1 T> 

0.8 - 



0.4 - \ 

0.2 - V P-- - 5 -..? 

o I . , . i , . , y ■ ■ ■ i ■ ■ ■ Y ■ ■ . i ■ ■ 

1 2 3 4 5 6 

Figure 3: Comparison of theoretical predictions (line) and experimental points for the ratio between 
the dressed and unperturbed precession frequency as a function of the amplitude of a square- wave 
like oscillating field. According to eq. ([3T|) the behaviour is described by the modulus of a sine 
function. 

A detailed analysis is performed in the case of a non-harmonic periodic transverse field 
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made of only two Fourier components (a fundamental one plus its k th harmonics). In this 
case the prediction (see eq. (|34p ) can be easily compared to the experimental results when 
varying one of the independent variables, as in the following. In Fig. H]we preliminarily show 
the behaviour of F as a function of the amplitude £ of the fundamental Fourier component 
and the relative amplitude r (see eq. (13T|) for the precise meaning of r, k and (p) of the k th 
harmonics, for fixed values of their relative phase and of k. 



k - 2 ip - 



k = 2 ip — 7r/3 





Figure 4: Graphs of the function F(£, r, k, ip), see eq. (|34p . for different values of k and (p. 

Due to technical limitations, only a few lobes of the ones shown in Fig. H] have been 
accessed to compare the model prediction with the experimental results. As an example, in 
Fig. the dependence of F on the amplitude ratio r is presented for several couples (k, £) 
of the frequency ratio and of the fundamental amplitude. 

One fact which stands out among the most interesting features of the system considered 
is that the non-linear nature of the observed phenomenon makes relevant the relative phase 
of the Fourier components of the oscillating field. Figs. and [7] are devoted to comparing 
the dependency of F on one of those parameters. For a large enough amplitude of both the 
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Figure 5: Comparison of theoretical predictions (line) and experimental points for different values 
of £ and A; as a function of r (see eq. (|3ip ). 



components, phase values exist causing the precession frequency to vanish. 

Two comparisons for k = 3 are made in Fig. El where differently located zeroes and local 
maxima of F appear. It is worth noting that, as confirmed by the experiment, in the case 
of k = 3 the dependence of the relative phase is symmetric with respect to the (p = tt value. 

Fig. [7] reports further comparisons made for k = 2. In these cases due to the small value 
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Figure 6: Comparison of theoretical predictions (line) and experimental points for k = 3, £ = 3 
and two values of and r for the ratio between the dressed and unperturbed precession frequencies 
as a function of the relative phase between the two Fourier components of the oscillating field (see 
eq. (J2U). 

of r, the interference-like pattern does not produce zeroes, while the number of local maxima 
is clearly dependent on the absolute amplitude £. For this value of k, both the model and 
the experiment show that the dependence of the relative phase has an additional symmetry 
with respect to the ip — 7r/2 value. 

The symmetries pointed out in the Figs. E] and can be interpreted in terms of properties 
of the Bessel function, according to the analysis presented in brief in the Appendix [A] 
As shown in that appendix, the symmetry around the 7r value is a consequence of a 2tt 
periodicity occurring for both odd and even values of k, while the additional symmetry 
around 7r/2 occurs only for the even values of k. 

The even-fc symmetry can be broken by the presence of a vectorial perturbation along 
the oscillation axis. In particular, a weak static magnetic field component added to the 
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Figure 7: Comparison of theoretical predictions (line) and experimental points for k = 2 and 
different values of £ and r for the ratio between the dressed and unperturbed precession frequencies 
as a function of the relative phase between the two Fourier components of the oscillating field (see 
eq. (J2U). 



alternating one results in an asymmetric curve 



w 



ip=0 



(38) 



ip=7T 



as can be seen in Fig. [SJ where it is also shown the theoretical simulations with and without 
the second order corrections (see eq. (126]) ). A look at the right panel reveals the importance of 
the corrections depending of the G n ,n^0 coefficients through the c quantity (see eq. (1251) ) 
in generating the asymmetry. 

This asymmetry has potential applications in designing and setting up apparatuses like 
the magnetometer used in our experiment. The scalar nature of our magnetometer makes it 
very sensitive to variations of the field component parallel to the static one (z component), 
while the transverse components can be detected and minimized with procedures responding 
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to the second order in £ - I n contrast, repeated measurements and equalization of Wq and 
w n constitute an effective procedure with a first order response in £ (see eq. (126]) ) for fine 
balancing the transverse dc component. In this perspective, an even-fc modulation technique 
may be proposed to achieve a vector response of this type of scalar sensors. 

It is also worth nothing that this asymmetry cannot be seen if the modulation is purely 
harmonic (see Appendix [B] for details), at least in the perturbation order considered in the 
described model. 
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Figure 8: (Upper) Comparison of theoretical predictions (line) and experimental points in presence 
of a small dc field superimposed to the ac one. The parameters are £o = yBx'/u) = 1.3 • 10 -2 , 
£ = 2,r = l,/c = 2 and ( = 0.14. The solid line is calculated using eq. (|26p while the dashed 
one using eq. (|22p . It is clearly visible the effect of a small £q- (Lower) The three second order 
corrections according to eq. ([26p . Notice how the term proportional to £q has the same value for 
(p = and (p = it and thus it is not responsible of the asymmetry. 
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V. CONCLUSIONS 



We have studied both theoretically and experimentally the effects produced on the Larmor 
precession by a transverse magnetic field, which oscillates at frequencies higher than the 
Larmor one, and has an arbitrarily large amplitude. The proposed model is suitable for the 
interpretation of the effect of a non-harmonically oscillating transverse field. In fact, the 
nonlinear nature of the system's response makes a harmonic expansion unfeasible and makes 
previously proposed analyses (considering only the harmonic transverse case) unsuitable for 
describing the general case of an arbitrary time-dependence of the transverse magnetic field. 
We have tested our model experimentally, studying in detail the case of a non-harmonic 
transverse field constituted by two Fourier components. The relevance of the relative phase 
of the two components in the determination of the observed dressed Larmor frequency was 
pointed out. We also demonstrated that the non-harmonicity of the alternating field gives 
rise to a first-order variation of the dressed Larmor frequency with respect to a small dc 
field superimposed to the alternating one. This observation is of interest for its potential in 
designing atomic magnetometers with vectorial response. 
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Appendix A: Properties of F. 

In this appendix we list some properties of the quantity (|34j) . An explicit expression can 
be obtained as 




(Al) 



where 



oo 




(A2) 



p=- 



oo 



20 



Using the property J_ m = (— l) m J m of the Bessel functions and some straightforward ma- 
nipulations one can see by inspection that 

Pk,s = (-l) (fc+1)s PM (A3) 

which in the case of odd k reduces to a trivial identity. For even k it is easily seen that 
Pk.2n+i = reducing the (p period in flAlj) from 2ir to tt. This feature is clearly visible in the 
Fig.0 



Appendix B: Asymmetry details 

Let's consider the second order corrections to Ql which give rise to the asymmetry. As 
can be seen from eq. f l26|) one needs to investigate the series ffMty . which can be rewritten as 

f~** I ' r* \ TUtJ \r~i |2 \r~i 12 



a 



Urg^Ur n , — (J- n ) ^ _ | <J"n \ ~~ -n\ (Bl) 



n * — ' n 

n=l n=l 



Then let's evaluate these expressions for the case of a general harmonic field as in eq. 
Notice that in real experimental situations the presence of the time lag tq is unavoidable. 
Thus we assume that tq is a random parameter uniformly distributed and at the end of the 
calculations all the quantities must be averaged. 

Using eq. (130]) and the property J_ n = (— l) n J n of Bessel functions one arrives at 



Re(a) = 2V ^^±1^ cos (2n + l)r b = (B2) 

whose average 



1 f 

Re(a) = — / Re(a) dr (B3) 
2tt J 

gives zero. We can thus conclude that a purely harmonic signal cannot give an asymmetric 
contribution. 

The non-harmonic case (see eq. (13T|) ) on the contrary shows the asymmetry. In fact let 
the general non-harmonic signal be 

0(t) = cos(r - t ) + rcos(A;(T - t ) + <p), (B4) 

using the same device of eq. ( 13"3"j) we can write 

G n (r Q ) = e*£( sin T 0+{r/k)sin(kT -ip)+(r/k) sin <p) Q -i ut (B5) 
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where the G n on the right-hand side can be read from eq. From this result easily 

follows that 

|G„(r )| 2 = |G n | 2 (B6) 

showing that the leading contribution to the Larmor frequency is unaffected by an overall 
phase. Moreover GL n cannot be simply related to G n (as in Bessel functions case), thus the 
b parameter in (1B1[) generally is non-zero. For the other series after some algebra one finds 

^f Re[G5(G B -G_ n )] Im[GS(G n -G-„)] . \ 

Re(a) = > < cosuto smnr > (B7) 

^— ' Ira n 

71=1 v y 

whose average is zero. 
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